On the theory of microwave absorption by the spin- 1/2 Heisenberg-Ising magnet 
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We analyze the problem of microwave absorption by the Heisenberg-Ising magnet in terms of shifted moments 
of the imaginary part of the dynamical susceptibility. When both, the Zeeman field and the wave vector of 
the incident microwave, are parallel to the anisotropy axis, the first four moments determine the shift of the 
resonance frequency and the line width in a situation where the frequency is varied for fixed Zeeman field. For 
the one-dimensional model we can calculate the moments exactly. This provides exact data for the resonance 
shift and the line width at arbitrary temperatures and magnetic fields. In current ESR experiments the Zeeman 
field is varied for fixed frequency. We show how in this situation the moments give perturbative results for the 
resonance shift and for the integrated intensity at small anisotropy as well as an explicit formula connecting the 
line width with the anisotropy parameter in the high-temperature limit. 



The total magnetization in an isotropic system of interacting 
spins rotates as a whole about the axis of a homogeneous ex- 
ternal field (see e.g. |jjy). We consider L spins-^, combining 

to a total spin S = S'^e^ + SyQy + S'^e^, S" = J^f^^ sf, 
in a magnetic field of strength h in z-direction. The Heisen- 
berg equation of motion for S with a Zeeman term —hS^ is 
solved by S(<) = {cos{ht)S^ +sin{ht)Sy)e^ - {sin{ht)S'= - 
cos{ht)Sy)ey + S^e^, a rotation counterclockwise about the 
z-axis. In ESR experiments this can be probed by circularly po- 
larized microwave radiation propagating along the z-direction. 
Since its wavelength is large compared to typical distances 
in spin systems, we may assume a magnetic field component 
of the form h(t) — A{cos{ujt)ex — sin(a;t)ej,), A > 0. It 
couples to the total spin as V{t) — —h°'{t)S" and produces a 
sharp resonance atoj ~ h (see Q below). If the spin system is 
perturbed by anisotropic interactions, this resonance is broad- 
ened, shifted or even split in a way that is characteristic of the 
microscopic interactions between the spins. 

For any spin system with Hamiltonian H linear response 
theory relates the observed absorbed intensity to the (imaginary 
part of the) dynamical susceptibility Hi] 

X+-{^,h) = — J^^dt&-\[S+it),S-])^. (1) 

Here = 5"^ ± 15'*' and (•)t stands for the canonical av- 
erage at temperature T calculated by means of the statistical 
operator p = Q-{H-hsn/T / ^^.^-iH-hs^/T _ Through this 
average the dynamical susceptibility depends on h and T. The 
absorbed intensity per spin, normalized by the intensity of 
the incident wave and averaged over a half-period tt/oj of the 
microwave field, is 

/(c.,/i) = |x+-(^,M- (2) 

In current ESR experiments in solids /(cj, h) is measured as 
a function of h for fixed lu. Of particular interest are experi- 
ments on quasi one-dimensional compounds (reviewed e.g. in 
in ID) which provide prototypical realizations of interacting 



many-body systems with strong quantum fluctuations. Still, 
the data are not always easy to interpret, because of a lack of 
reliable theoretical predictions. 

Most of the existing theories are based on a priori assump- 
tions about the line shape and typically apply for limited ranges 
of temperature and magnetic field. Field theoretical approaches 
[1 1 are restricted to small temperatures and small (but not too 
small) magnetic fields. The more traditional approaches IJUS] 
rely on the high temperature approximation. Purely numerical 
approaches |6| are unbiased, but the extrapolation of the data 
to the thermodynamic limit of large chains may be difficult. 

A remarkable result for the resonance shift in one- 
dimensional antiferromagnetic chains, being valid at arbitrary 
temperatures, was obtained in (|2|. It utilizes the exact nearest- 
neighbor correlation functions of the isotropic spin-^ Heisen- 
berg chain. Here we present an alternative framework for the 
derivation of the resonance shift which, in the limit of small 
anisotropy, reproduces |7|. In our approach the anisotropy is 
treated non-perturbatively, and it allows us to derive an exact 
formula for the line width 'in frequency direction' at fixed 
magnetic field, as well as a new explicit expression for the 
ESR-line width in the high temperature regime. 

We consider an important example of anisotropic interac- 
tions described by the Heisenberg-Ising Hamiltonian 

H = jY^is^s^ + s\s) + (1 + 5)sls^) . (3) 

Here the sum is over nearest neighbors, and 5 is the anisotropy 
parameter. If (5 = 0, then ([T]l, Q imply that the normalized 
absorbed intensity is 

/(w, h) = tt5{uj ~ h)hm{T, h) . (4) 

It is proportional to the magnetic energy hm{T, h) per lattice 
site. This case includes the familiar paramagnetic resonance 
(Zeeman effect) for J = 0, for which the magnetization is 
m(T, h) = ^tl^i^). 

For non-zero 6 the function _ is unknown and hard to 
calculate. Still, some more elementary spectral characteristics. 
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FIG. 1. Resonance shift Suj/ J for the ID model in the critical regime 
at 5 = —0.1 as function of the magnetic field. Crosses from fully 
numerical calculation for a finite chain Hamiltonian of 16 sites. 
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FIG. 2. Line width Aco/J for the ID model in the critical regime at 
(5 = —0.1 as function of temperature. Crosses from fully numerical 
calculation for finite chain Hamiltonians of 16 and 24 sites. 



such as the position of the resonance or the line width, may be 
expressed in terms of certain static correlation functions that 
determine the moments of the normalized intensity function. 

Let us assume for a while that our chain is large but finite. 
Then the spectrum is bounded and the integrals 



In = 



(5) 



exist for all non-negative integers n. Since h) is non- 
negative everywhere and since /q > 0, we may interpret 
/(w, h) //o as a probability distribution and the /„ as its mo- 
ments. As we shall see, it is convenient to express the /„ in 
terms of another closely related sequence of integrals 



TO„(T, h) = J 



2ti 



ioj-hrx'i-i^.h) (6) 



which, by slight abuse of language, will be called (shifted) 
moments. Again they exist for every finite chain. 
By definition the shift of the resonance for fixed h is 



6lj 



h 
lo 



J 



Jm2 + hmi 
Jrrii + hniQ 



(7) 



A measure for the line width is the mean square deviation 



Iq J mi 



hniQ 



So:' 



(8) 



Hence, in order to calculate the resonance shift and the line 
width, we need to know the first four shifted moments mo, nii, 
TO2, TO3 of the dynamic susceptibility X+-- 

In the following we shall employ the notation adx ■ = [X, ■] 
for the adjoint action of an operator X. Then = 
^-iht^it^dHS+^ since [H, S''] = and [S'^ S+] = S+, and it 
follows with ([T]i and (|6| that 



,atlJ„S-|> 



T ■ 



(9) 



The latter formula shows that the moments m„ are static 
correlation functions whose complexity grows with growing n. 
The first few of them can be easily calculated. We shall show 
the results for the one-dimensional model, for which 



mo^^{[S+,S-])^ = ^{S')^, (10) 



which is the magnetization per lattice site. The subsequent 
moments are less intuitive. 



nil = ^{•''t^2 ~ 2siS2)t , 



(11a) 



m2 = ^S^{sl+'^444- '^44s3)t , (lib) 



sjs^ + 8sls2S^sl + 2slsl - 4:sls2 
+ 2s+S2--8sfs^))^ 



(11c) 



They are certain combinations of static short-range correlation 
functions. This implies, in particular, that they all exist in the 
thermodynamic limit L — > 00. It follows that the line shape 
for fixed h cannot be strictly Lorentzian, as is often assumed 
in the literature. 

All static correlation functions of the one-dimensional 
Heisenberg-Ising model are polynomials in the derivatives 
of three functions u, lo' and ip fS) which, as is common in 
integrable models, can be expressed in terms of the solutions 
of certain well behaved linear and non-linear integral equa- 
tions ||9|- This is the reason why in this case the moments 
toqi JTii, r7i2, ?7i3 can be calculated exactly by means of the 
techniques developed in Ll^ . 

We parameterize the anisotropy as 5 = ((7 — l)-^/2g. Then, 
with the shorthand notations (y9(„) = d''^tf[x)\x=Q, f{m.n) — 
d^d^f{x,y)\^^y^o^ for / = uj,uj', we obtain 
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1 _ (<? - +4q+ 1)^1^^^^ (^3 _ 



'(0,0) 



mo = --^(o), mi ^g^, ^^^^^^^ , 

(q - 1)2 

^2 = ^^g^4 [4g(g + l)(g^ - l)(w(o,2)</'(o) - 2w(i^i)(^(o) - a;(o,o)'y3(2)) 

+ (g2 - 1)2(^2 ^ 4^ ^ l)(^^_2)</'(0) + '^(o.i)¥'(2)) - 16g'(g - l)V(o)] , 

'^^ 98304q8(g4 _ i)(^6 _ 1) 

[iGq^q^ ~ l)3(g^ - 1)(96 - l)iq^ + 4g + l)(t^(o,2)C^(o,i) + ^(0,0)^(1,2)) 

+ 64g2(q2 _ 1)4(2^10 _ ^9 ^ 4^8 _ _ ^3^6 - iV - 129"^ - Aq^ + Aq^ - q + 2)a;(o,o)^^(i,i) 

- 16g2(q2 _ 1)2 (g4 _ ^)(^6 _ i)(3^4 ^ 14^2 ^ 3)^,^^^^ 

+ 8q\q' l)(g4 - Ifiq' -l){q + l)^(8^[i,2) " ^(2,3)) 

+ 192g''(92 - 1)2(^4 _ -^)(^6 _^ ^8^4 _^ 8^3 _^ ^8^2 ^ 1^^ ^) 

+ 64g2(q2 _ i)2('g4 _ -^)(^ ^ l)^(2g* - + 26^^ - 49g^ + 28^"^ - 49^^ + 26q'^ - 5q + 2)w(i,i) 

- 16g2(q2 _ i)2(g4 _ ^)(^ ^ ^)2(^8 _ ^7 ^ ^6 ^ ^5 + 2q^ + q^ + q^ - q + l)(2c^(i,3) - 3w(2,2)) 
+ 64g2(q4 _ i)(^6 _ -^^(3^8 ^ ^ 24g5 _ i30g4 + 24q3 + 2g2 + 3)tj[o_i) 

+ (g^ - l){q^ - l){q + lf{q^" - 2q^ + 25q^ + 16g^ + 1189*^ + 164?^ 

+ 118g* + 16g3 + 25g2 - 2q + l)(w[o,3)W(i.2) + '^(0,1)^(2,3)) 

- 15369^(5-^ - l)(4g« - - 2q^ - 6q^ + 8q^ - 6q^ ~ 2q^ - 9g + 4)a;(o,o) 

+ 4^2(^6 _ + 1)2 (q2 _^ - 2g^ + 32^*^ + 50q^ + IQq^ + SOq^ + 32^2 - 2q + 5) 

(2^(1,3)^(0,1) ^ 3i:J(2,2)W(o,i) +^(0,2)1^(0,3) ^ 2w(i,i)W(o,3) ~ 3W(o,2)W(i_2) ~ "^(0,0)^(2,3)) 

- 16g2(q + \f{q^'^ - q^^ + 'iq^^ + 99^^ + 47ql2 ^ 45^11 ^ gggio ^_ gi^S _^ 128^8 _^ gi^7 ^ gg^e 

+ 45g^ + 47g'' + 9g^ + 8g^ - 9 + l)(3a;2g 2) - 6w(i^i)W(o,2) + 2^(0,0)^^(1,3) - 3w(o,o)W(2,2))] ■ (12) 



These functions represent the moments in the thermodynamic 
Umit. Since they can be calculated to arbitrary precision, we 
obtain numerically accurate results for the resonance shift and 
for the line width as functions of temperature or magnetic field 
over the whole range of the phase diagram. In particular, our 
approach is not restricted to small anisotropics. Examples for 
8 = —0.1 are shown in figures[T and|2] We find a broadening 
of the line width as defined by ( 8]) for small temperatures in 
the critical (S < 0) as well as in the massive (6 > 0) regime 
(latter case not shown here). 

At first sight this seems to contradict experimental results 
m which claim a narrowing. Still, one has to take into ac- 
count that usually in the analysis of experimental data rather 
different definitions of the line width, as e.g. the distance be- 
tween the tm-ning points right and left to the maximum of the 
intensity, are used. In particular, if the intensity distribution 
has long shallow tails the definition ([8]l will give considerably 
larger values than the distance between the turning points. We 
have performed a detailed numerical study of the dynamical 
susceptibility for chains of 16, 20 and 24 lattice sites (to be 
published elsewhere) and we see indeed such long tails at low 
temperature (compare also [6|). In experiments they may be 



misinterpreted as background stemming from couplings of the 
spin chain to other degrees of freedom, but in fact they are due 
to the spin-spin interactions and are part of the true absorption 
line. For the determination of the resonance shift tails are ex- 
pected to have less influence. As long as they are not too much 
asymmetric the shift calculated by means of (|7| should agree 
with the shift of the maximum of the absorbed intensity. 

In current ESR experiments the microwave frequency lu is 
kept fixed and the Zeeman field h is modulated. This means 
that, as opposed to most of the theoretical treatments, including 
our considerations above, the absorbed intensity I{oj, h) ~ 
wx"_(w, h)/2 is determined as a function of h for fixed 
and the resonance shift and line width are measured in '/i- 
direction'. Away from the isotropic point {5 = 0), where 
h) is symmetric and the absorption line is extremely 
narrow, this may clearly lead to rather different values. To be 
closer to present-day ESR experiments one should calculate 
resonance shift and line width in terms of the moments of the 
dynamical susceptibility in '^.-direction'. 

We define 



M„(T,c^) - J- 



dh 
2n 



(13) 
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For these functions we obtain the representation 

NU{T,u) = (-1)"^ ^4^^jr,..) , (14) 

where the superscript (fc) denotes the /cth derivative with re- 
spect to the second argument. We see that the Mn are deter- 
mined by infinitely many of the m„ and their derivatives, i.e., 
they depend on static correlation functions for arbitrarily large 
distances. For this reason they cannot be calculated by our 
exact method above. Yet, in certain cases only finitely many 
terms of the series are needed for a good approximation. 

We first of all express the resonance shift 5h = (h) — uj and 
the mean square deviation from the center of the absorption 
peak A/i^ = (/i^) - (h)"^ in terms of the M„, 



with 



6h 



Ml 
Mo 



M2 
Mo 



Ml 
Mi 



(15) 



We have identified two cases, where these formulae simplify 
and finitely many of the m„ are enough to determine Sh and 
Ah approximately. 

The equation for the resonance shift simplifies for small 
anisotropy \6\ <^ 1. Since AIq — itiq + 0{S), Mi — mi + 
0{5^) and, generically, nii itself is of order 6 (see ( 1 1 1, ( 14 1) 
we obtain to linear order in 5 



6h 

T 



mi 
mo 



(16) 



In ISl IT] the same equation was obtained by a more intuitive 
reasoning. It leads to results which compare rather well to 
experiments |7|. However, some care is necessary with the 
interpretation of ( 16 1, rni/6 vanishes at 5 = /i = 0. It follows 
that mi = 6{ah+b5+. . . ) with some coefficients a, b, whence 



h must be large compared to (5 for ( 16 1 to be applicable. Note 
that all higher moments m„ are of order 6^. Hence, there 
is no simplification for small anisotropy, like in ( |16| l, for the 
line width. But the integrated intensity Mo has again a finite 
approximation to first order in S, Mq = mo — Jm'i . 



The representation ( 14 1 is a series in ascending powers of 



J/T (with still temperature dependent coefficients). This can 



be used to evaluate ( 15 1 asymptotically for high temperatures 
It turns out that that the leading terms in the J/T expansion 
of mi and Jm'2 cancel each other: Sh 
high-temperature limit T ^ J, and 



^6 in the 



Ah 



V2 



(17) 



for arbitrary microwave frequency uj. This formula provides a 
simple means to directly measure the anisotropy parameter 5. 

It may be instructive to illustrate our formula with one of the 
few explicit results, namely with the formula for the intensity 
in the free fermion case 5 = — 1 at T oo 1 11 1. In this case 
I{uj,h) ^ (w^/ J^) exp(— (cj — h)^/J^) and, in agreement 



0, 



we obtain the line width ^ 



whereas the 

V 2 



width m omega direction depends on h. 

Our work is the first exact result for the resonance shift 
and the line width in microwave absorption experiments on 
the Heisenberg-Ising chain. The reduction to moments is not 
restricted to the integrable case and may be interesting for the 
two- and three-dimensional models as well. Our approach is 
unbiased. It makes no a priori assumptions about the shape of 
the spectral line. As opposed to all other approaches it is valid 
for all temperatures and magnetic fields and in addition for 
arbitrary values of S. For small 5 close to the isotropic point 
we recover the result of [71 for the line shift. The resonance 
shift 5uj / J ox 5h / J and the line width Alu / J ox Ah/ J defined 
in terms of moments show a simple scaling behavior. They 
depend on the exchange interaction only through the ratios 
T / J and h/ J. In this sense the curves in figures [T] and |2] are 
universal. 

The intensity /(cj, h) is a function of uj and h. With our 
definitions of Suj and Ao; we determine the resonance shift and 
the line width in w-direction as functions of h, while in standard 
ESR experiments the resonance shift 5h and line width Ah in 
/i-direction are measured as functions of lo, which should be 
clearly distinguished. For the resonance shift it follows from 
Q, ([T6]l that 5h{T, u) = -Slu(T, h)\h=cj to hnear order in S. 
For the line width there is no such simple relation between 
Au! and Ah. However, for Ah we obtained the simple high- 



temperature formula ( 17 1 which we suggest to be useful to 
measure the anisotropy directly. We are further convinced that 
it may be worth trying to measure Au, which is now known 
exactly, directly in multi-frequency ESR experiments. 
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